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ABSTRACT 

In this paper we discuss the evolution of gravitationally unstable pre-galactic discs that 
result from the collapse of haloes at high redshift z « 10 or so, which have not yet 
been enriched by metals. In cases where molecular hydrogen formation is suppressed 
the discs are maintained at a temperature of a few thousand degrees Kelvin. However, 
when molecular hydrogen is present cooling can proceed down to a few hundred degrees 
Kelvin. Analogous to the case of the larger scale proto-galactic discs, we assume that 
the evolution of these discs is mainly driven by angular momentum redistribution 
induced by the development of gravitational instabilities in the disc. We also properly 
take into account the possibility of disc fragmentation. We thus show that this simple 
model naturally predicts the formation of supermassive black holes in the nuclei of 
such discs and provides a robust determination of their mass distribution as a function 
of halo properties. We estimate that roughly 5% of discs resulting from the collapse of 
haloes with M w 10 7 M Q should host a massive black hole with a mass Mbh ~ 1O 5 M . 
We confirm our arguments with time-dependent calculations of the evolution of the 
surface density and of the accretion rate in these primordial discs. The luminosity of 
the outer, colder disc is expected to be very low (in the range of a few thousand Lq), 
while the formation of the black hole is expected to produce a burst with a luminosity 
of a few times 10 9 L Q . This mechanism offers an efficient way to form seed black holes 
at high redshift. The predicted masses for our black hole seeds enable the comfortable 
assembly of 10 9 M Q black holes powering the luminous quasars detected by the Sloan 
Digital Sky Survey at z = 6 for a concordance cosmology. 

Key words: accretion, accretion discs - black hole physics - galaxies: formation - 
cosmology: theory - instabilities - hydrodynamics 



1 INTRODUCTION 

The local demography of black holes at the centers of galax- 
ies suggests that bl ack hole formation is a generic feature 
of galaxy formation iMag orrian et aJll998tlTremaine etafl 
l2002l:lFerrarese fc Merrittl2000l) . Mergers and accretion pro- 
cesses, in fact, merger induc ed accretion are implicated in 
the assembly of black holes ( Kauffmann & Haehnelt 2000; 
IVolonteri etaD l2003t iDi Matteo et alll2005l) . However, all 
models require the formation of seed black holes at high 
redshift. The inferred large masses (M bh ~ 1O 9 M ) of 
the black holes powering luminous quasars at z ~ 6 de- 
tected by the Sloan Digital Sky Survey (SDSS) offers a 
challenge to mechan isms for product ion of seed black holes 
at higher redshifts <Fan et al.ll2004f) . In a ACDM cosmol- 
ogy, these black holes have to assemble such large masses 
within a Gyr. The popular picture is that the remnants 



of the first generation of metal-free stars that form with 
an initi al mass function that is biased to wards higher 
masses jAbel et alJ l200d : iBromm et ail 12002?) provide the 
initial seed bla ck holes with masses of the order of a 
IP's to 100 jMadau fc Reedl200lt I Volonteri et alJl2003t 
iRicotti fc Ostrikerl2004l:lMapelli et all200dh Growing these 
seeds to the requisite masses by z — 6 to match the abun- 
dance of observed SDSS quasars requires various kinds of 
fine-tuning including phases of maximal grow t h pow ered 
by super-Eddington accretion. IVolonteri fc R ccs (2005) con- 
sider the growth of black holes in metal-free halos with 
T > 10 4 K that cool via atomic hydrogen lines leading 
to the formation of a fat, gas disc. They argue that an 
episode of super-Eddington accretion assembles black holes 
with masses of 10 6 M at z ~ 15 - 20. 

There are alternative models that predict the direct 
formation of more massive seeds with masses of about 
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10 5 M Q . These ranee from scenarios based on the formation 
of supermassive objects formed dir e ctly o u t of the collapse 
of d en se gas (|Haeh nelt_fcFieesj [1233; lUmcmura et_alj 
19931: lLoeb fc Rasiol Il994t lEisenstein fe Loebl 1 1993 
Bromm fc Loet 12001 iKoushiappas et alJ booi 
Begelman et al, l200fift . The key limiting factor for 
these models is the dispo sal of the angular mo men- 
tum. One approach taken bv lEisenstein fc Loebl lll995T) and 
IKoushiappas et alJ <l2004ft is to argue that it is preferen- 
tially either low spin halos with consequently low angular 
momentum gas or the low spin tail of the gas distribution 
in halos that can cool efficiently that form these seeds. 
Significant transfer of angular momentum is still required 
to go from the central massive object to the final collapsed 
hole that is expected to occur via the post-Newtonian 
instability. No n-axisymmetric struc tures like bars have 
been proposed iShlosman et alJll99d) to efficiently transfer 
angul ar momentum at these late stages. In a recent paper, 
iBeeel man et all ]2006ft have pursued this picture and argue 
that a low specific entropy 'quasi-star' is produced as a 
result of the cascade of the bars-within-bars instability 
leading to the formation of a black hole of a few solar 
masses. 

In this paper, we argue for the formation of more mas- 
sive black hole seeds (~ 10 j A/q) from gas cooling in pri- 
mordial halos via the growth of gravitational instabilities. 
We offer a picture wherein accretion processes and fragmen- 
tation criteria are addressed in a coupled fashion. The out- 
line of our paper is as follows: in Section 2, we outline the 
scenario for pre-galactic disc formation in dark matter ha- 
los and discuss their accretion properties, the fragmentation 
criteria which determine the fate of the gas are discussed in 
Section 3 and the implied mass distribution of central con- 
centrations is derived in Section 4. Detailed time-dependent 
models that confirm the analytic calculations in the previous 
sections are presented in Section 5, the resultant luminosity 
and potential observability of these sources is discussed in 
Section 6 followed by a discussion and implications of our 
work. 



2 PRE-GALACTIC DISC FORMATION 

In the context of the hierarchical, cold dark matter dom- 
inated structure formation paradigm, the baryonic compo- 
nents of galaxies assemble upon dissipation within dark mat- 
ter halos. Generic to these models is the assembly of discs 
in these haloe s as a consequence of dissipational co llapse 
of the baryons (White & Reesjl978tlFall fc Efstathiolll980l: 
iMo et al.ll998t[Dalcanton et alll997ft . Here, we revisit some 
of the standard arguments put forward to describe the for- 
mation of proto-galactic discs inside a spherical dark matter 
halo of mass M. In partic ular, we adopt the formalism devel- 
oped bv lMo et alJ <ll998ft . and extend it to include the effect 
of angular momentum redistribution and mass accretion due 
to gravitational instabilities in the gas disc. 

The arguments that form the basis of the standard disc 
formation picture are simple and elegant. It is assumed that 
a fraction md of the total mass of the halo, containing a 
fraction ja of the angular momentum of the halo (usually 
assumed to be equal to md), collapses and forms a disc. 
Note that the baryonic mass in the halo is assumed to be a 



fixed fraction of the halo mass, and the proportionality fac- 
tor is taken to be the universal baryon fraction /(,. Assuming 
a radial profile for the resultant gas density distribution in 
the disc (for example, exponential) enables the calculation 
of its key physical properties: the disc scale length, central 
surface mass density, and crucially its stability under self- 
gravity. Star formation is then assumed to occur for those 
discs that are gravitationally unstable. The results of ear- 
lier work based on this paradigm de pend on the det ailed 
assumptions made. For example, while lMo et al] lll998ft con- 
sider both an isothermal sphere and a NFW ^Navarro et al] 
11997ft profile for the dark mat ter halo and a globa l criterion 
for gravit ational instability. | Oh fc Haimanl ]2002ft consider 
an NFW ]Navarro et aljll997ft halo and a local instability 
criterion, based on the Toomre-Q parameter. These treat- 
ments present 'static' models for disc assembly. In this work, 
we examine the fate of these pre-galactic discs taking their 
time dependent evolution into account. 

The limitation of previous analyses is that they consider 
the properties of the discs a posteriori, after all of their mass 
has been assembled, and do not consider the evolution of the 
discs during their formation process. They have argued that 
for some values of the relevant halo parameters, in particu- 
lar, for low values of the spin parameter of the halo A 1 , the 
disc can be self-gravitating. In principle, these analyses also 
allow the disc to have very low values of Q or equivalently to 
be strongly unstable in terms of a stability criterion. How- 
ever, it is important to note here that as soon as the disc 
becomes massive enough to be marginally stable, it will de- 
velop structures that will re-distribute angular momentum 
and mass through the disc, preventing the surface density 
from becoming too large and the disc from becoming too un- 
stable. A suggestion that this kind of matter redistribution 
in primordial discs can lea d to the formation of s eed black 
holes was ma de recently b y Begelman^talJ <l200a) and also, 
previously, bv IKoushiappas et al] ]2004ft but using substan- 
tially different arguments compared to our work here. In 
practice, the accretion of matter to the center will stop as 
soon as the surface density becomes low enough to make the 
disc stable, so that in the final state the disc will be exactly 
marginally stable. We explore this scenario and the conse- 
quences of angular momentum redistribution in detail here. 
In this section, we propose some simple arguments that al- 
low us to estimate the final disc mass and the mass that will 
flow into the central region, leading to the formation and 
growth of a seed black hole. In Section |5] we support these 
simple arguments with full time-dependent calculations. We 
are therefore able to keep a detailed, self-consistent inven- 
tory of the fate of the gas. 

We assume that in the final configuration, once the pre- 
galactic disc has assembled in the dark matter halo, only a 
fraction rrtf of the total halo mass is retained in the disc 
(note that we track the baryons that remain in the disc and 
those that are accreted towards the center separately), while 
the remainder m a = m-d — m,f is accreted towards the center 
of the galaxy providing fuel for a growing black hole. The 
final angular momentum per unit mass of the disc, however, 



The spin parameter A of a dark matter halo is defined as A = 
J \E\ 1 /' 2 /G M 5 / 2 , where J is the angular momentum of the halo, 
E its total energy and M the halo mass. 
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remains equal to the initial value j d , since, in any accretion 
process, angular momentum is transported outwards, and 
the only loss of angular momentum due to accretion is re- 
lated to the advection of it out of the inner radius of the 
disc. For an isothermal dark matter halo p(r) oc 1/r 2 with a 
flat rotation curve, the angular momentum advected inwards 
due to the accretion of a mass m a is j a = jdRin/Rd << jd 
(here Ri n and Rd are the inner and outer radius of the disc, 
respectively). It is assumed here that the disc is embedded 
in a spherical, isothermal, dark matter halo, with virial tem- 
perature T v i r and a constant rotational velocity Vh- To eval- 
uate the stability of the disc, we use the Toomre parameter 
Q which is defined as: 
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where c s is the sound speed, k — y2Vh/R is the epicyclic 
frequency, R is the cylindrical radial coordinate, and E is the 
surface mass density of the disc. We consider here the earli- 
est generations of discs, which have not been metal-enriched 
and therefore are able to cool only through hydrogen. We 
present in the Appendix the detailed structure of these discs. 
Here we note that, in thermal equilibrium, if molecular hy- 
drogen formation is inhibited, these discs are expected to be 
nearly isothermal at a temperature of a few thousand K. 
If molecular hydrogen is present, however, the disc can be- 
come much colder, with a temperature of a few hundred K. 
A marginally stable, isotherma l disc has th e same surface 
density profile as a Mestel disc jMestellll963l) . We therefore 
assume that: 



(2) 



Dark matter halos are characterised by their mass M and 
their angular momentum J, or equivalently by their spin 
parameter A. The distribution of spin parameters for cold 
dark matter halos as determined from cosmological N-body 
simulations is given by: 
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where the mean spin A = 0.05 and the dispersion is ax = 0.5 
JWarren et alll992l) . Requiring that the mass of the gas disc 
is a fraction (nid — m a ) of the halo mass and that its angular 
mome ntum is a fraction jd of the halo, following iMo et all 
(1998) we easily get: 
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where V2oo = Vh/WH(z) is the halo virial radius (the radius 
within which the density is 200 times the critical density of 
the Universe) and H(z) is the Hubble constant as a func- 
tion of redshift z. We can now evaluate Q, given the disc 
temperature T gas : 
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We use the fact that the final disc configuration is exactly 
marginally stable, and thus obtain m a by requiring that Q = 



Q c , a critical value above which the disc is gravitationally 
stable and no accretion can take place, due to the lack of an 
angular momentum transport mechanism. This gives: 
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Note that, although the result summarized in eqn. Q 
is based on the assumption that angular momentum is redis- 
tributed within the disc, it does not depend on the specific 
viscosity mechanism, and is therefore very robust. The only 
assumption made is that whatever this mechanism is, it is 
only active when the disc is gravitationally unstable (a very 
reasonable assumption for such primordial discs, for which 
the main source of viscosity comes from gravitational insta- 
bilities). As shown, the amount of mass that will be con- 
centrated in the central regions of these pre-galactic discs 
depends only on halo properties (such as the spin parame- 
ter A and the fraction of baryonic mass that collapses to the 
disc m d ), the ratio between gas temperature and halo virial 
temperature, and on the threshold value of Q, which has a 
very small range of variation around Q c w 1. 

Equation J7J then provides a powerful link between the 
properties of dark matter haloes and the mass of massive, 
seed black holes that can grow within them. This implies 
a larger mass concentration in the center for haloes with 
low spin parameter A or with high virial temperature T V1I 
(corresponding to higher halo mass mass M). The fact that 
seed black holes are likely to for m preferentially in low spin 
halos w as also demonstrat e d bv lEisenstein fc Loebl dl995f) 
and by iKoushiappas et alJ i2004f) . using a different set of 
arguments. However, for high mass haloes, there is the added 
possibility that the disc fragments. We discuss the issue of 
fragmentation in Section|H]and we postpone to Section|l|the 
discussion of the mass distribution of central concentrations 
given the distribution of the spin parameter A. 



3 FRAGMENTATION CRITERIA 

The arguments outlined in the previous section implicitly 
assume that, when the disc is gravitationally unstable, it is 
able to redistribute mass and angular momentum so as to 
accommodate the incoming mass inflow from the halo, main- 
taining a marginally stable state. Actually, in some cases, 
this is not going to be the case and the disc might undergo 
fragmentation. The criterion for the fragmentation of a grav- 
itationally unstable disc has been studied extensively in the 
past few years (especially in relation to models of planet for- 
mation) . lOammieT l)200 jl has shown that being unstable (i.e. 
having Q < 1) is a necessary but not sufficient requirement 
for disc fragmentation. He has shown that fragmentation 
occurs only when the cooling time is faster than the disc 
dynamical timescale. An equivalent fragmentation criterion, 
that directly relates ac cretion and fragmentation properties, 
has been discussed bv lRice et, alJ l|2005h . Indeed, they have 
shown that the gravitationally induced stress in a thin disc 
cannot exceed a critical value, which, measured in terms of 
the standard a description of transport properties of accre- 
tion disc, corresponds to a c w 0.06. If the disc is required 
to provide a larger stress, by either strong cooling or by a 
large mass inflow, it will undergo fragmentation within a dy- 
namical time-scale. In the context of pre-galactic discs that 
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we are considering here, this will then lead to efficient star 
formation in the disc. 

In the case considered here, the constraints imposed by 
cooling are not important, since, due to the steepness of the 
cooling function when the gas is predominantly hydrogen (in 
either the molecular or atomic form), slight changes in the 
temperature can easily bring the disc into thermal equilib- 
rium. The constraint coming from the inflow of mass from 
the halo is much more important. The mass accretion rate 
Mh from the halo can be simply estimated by requiring that 
a mass m d M (where M = V^r2oo/G is the total halo mass) 
collapses in a free-fall time-scale tg = rzoo/Vh- We therefore 
get: 



M h = m d 



(8) 



For a self-gravitating disc the standard mass conserva- 
tion equation for a steady disc: 



2ttz/E = M, 



(9) 



where M is the mass accretion rate and v is the viscos- 
ity (described here with the standard a-parameterization) , 
takes the simple form: 



M = 2a- 
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from which we see that in a steady state any isothermal disc 
will be characterized by a constant viscosity coefficient a. 
Note that in eqn. there is a factor 2 rather than the 
standard factor 3 of Keplerian discs, because the rotation 
curve here is assumed to be flat. 

What would happen to the transport properties in the 
disc once the fragmentation threshold is overcome is still 
not clear. A plausible assumption is that the accretion rate 
would be capped at its threshold value, with extra heating 
needed to reach thermal balance. This heat could be pro- 
vided by the star formation process. The disc would then 
not be able to transfer down to its inner parts an M larger 
than: 



M max = 2a c -^, 

G 



(11) 



where a c « 0.06. 

If we consider the case of cooling via atomic hydro- 
gen and assume that T gas m 4000 K, we find that M max ~ 
lO _2 M0/yr. In contrast, for the cold discs arising when 
molecular cooling is efficient, the maximum accretion rate 



is significantly lowered, down to M n 



3 x lO _4 M /yr 



If the disc evolves to a steady state, one should expect 
that the mass flow in the disc would equal the inflow rate 
Mh. Actually, as the disc builds up, only a fraction of the 
matter flows into the inner disc, since the remaining has to 
be transported outwards to take up the angular momentum 
lost by in-flowing material. Since the viscous time-scale in 
the outer disc is very long, even if the inner disc is almost 
in a steady state, the inner flux Mi n will still be a fraction 
of Mh, that we will calculate explicitly in Section We 
anticipate here the result that the mass flux in the inner 
disc is at most a fraction / = 1/2(1 + m a /md) of the total 
inflow of matter. Note that 1/2 < / < 1, so even in the 
extreme case where no mass in accreted (i.e. m a = 0), the 
mass flow is only reduced by a factor of a half. We expect 
fragmentation if: 



M in = m d ^i + >2a c ^. (12) 

G 2 \ m d / G 

We can rearrange eqn. I12H in terms of the ratio of the virial 
temperature to the gas temperature: 
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Assuming a c = 0.06 and m d = 0.05, we can therefore iden- 
tify three possible behaviours: 

• Tvir/Tgaa > 2.9 : these haloes will form discs that will 
fragment and form stars. 

• 1.8 < TVir/Tgas < 2.9: these haloes will form fragment- 
ing discs only if the spin parameter is significantly low, so 
that m a /m d ~ 0{\). 

• T v i r /T gas < 1.8: these haloes will not produce frag- 
menting discs, but since the halo mass is relatively low, the 
amount of mass that can be concentrated in the center will 
also be correspondingly low (cf. eqn.|7|l. 

In practice, whenever it happens, fragmentation will act to 
suppress the concentration of the largest amount of mass 
in the center, whereas haloes with relatively smaller central 
concentrations are by and large unaffected by fragmentation. 

The criterion for fragmentation described above, and 
the critical value of a, have been obtained in the context 
of thin and Keplerian discs. The haloes that give rise to 
large concentration of mass in the center have T v i r > T gaa 
and therefore the discs that form within them are relatively 
thick. In this case, the Jeans mass in a fragmenting disc 
is only slightly smaller than the total disc mass, so that 
fragmentation might be less likely in this case than for thin 
discs, effectively leading to an increase of the critical value 
of a. 



4 MASS DISTRIBUTION OF CENTRAL 
CONCENTRATIONS 

We now discuss the distribution of central mass concen- 
trations, given a distribution of halo spin parameters. We 
determine the distribution of the mass available to feed a 
growing black hole in the center of the pre-galactic disc, 
Mbh = m. a M, based on eqn 0, while the distribution of A 
is given by eqn. @. We also need to specify the values of m d 
and Q c . We assume m d = 0.05, but we also consid er the case 
were m d = 0.1. As customary, we take j d = m d iMo et alJ 
1998). The critical value for the Q-parameter, Q c is more 
uncertain. A self-gravitating disc is marginally stable to lo- 
cal, axi-symmetric instabilities at Q = 1. However, non axi- 
symmetric and global instabilities cause the disc to become 
unstable at larger values of Q. The prescription described 
above implies that at Q = Q c the disc provides no stress, 
while on the other hand, lLodato fc Ricd i2004l 120051) have 
shown that a Q = 1 disc is able to deliver a sizable stress. 
We expect Q c to be of the order of, but relatively larger 
than, unity and therefore consider two cases: (i) Q c = 2 and 
(ii) Qc = 3. 

4.1 The case of atomic hydrogen cooling 

First, we consider the case when molecular hydrogen for- 
mation is inhibited and the gas can consequently cool down 
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Figure 1. Cumulative distribution of the mass available to the formation of a seed black hole: The parent haloes are at a redshift z = 10 
and we assume here that m d = 0.05. The line styles refer to M = 3 X 10 7 A/ Q (solid), M = 2 X 10 7 A/ Q (dotted), M = 1O 7 M (dashed). 
The dot-dashed line refers to the case of an M = 3 X lO 7 A/0 halo, when fragmentation at the lowest A is allowed. For the two lower halo 
masses fragmentation will not occur. 



to T gas « 4000 K through atomic hydrogen cooling only 
(see Appendix for details). As discussed above, for a given 
r gas , haloes with larger T v j r produce larger central mass 
concentrations. However, we also have to require that the 
disc does not fragment, which then limits the analysis to 
Tvir ^ 3T gas « 10 4 K. The typical redshift of collapse of 
haloes with this virial temperature is around z « 10. 

In Fig. Q we show the results for the two cases Q c — 2 
and Q c = 3, respectively. For both cases, the haloes are 
assumed to be at a redshift z = 10 with mass in the range 
of 1O 7 M0. The lines refer to 3 values of the halo mass, 
M = 3 x 10 7 M Q (solid line), M = 2 x 1O 7 M (dotted line) 
and M = 1O 7 M (dashed line). While for the two lower 
mass cases T v i r < 1.8T gas and according to the results of 
the previous section, fragmentation is not expected to take 
place, for the highest mass case the disc might fragment for 
the lowest values of A. The dot-dashed line refers to the case 
where the possibility of fragmentation is included and in this 
instance, we assume that these discs simply do not produce 
any central mass concentration. 

It can be seen that, especially for the Q c = 3 case that 
a significant fraction of discs host central mass concentra- 
tions. For dark matter haloes of mass M = 3 x 1O 7 M , 
approximately 12% of discs have a central concentration of 
mass larger than 10 M , while the same mass concentration 
can be achieved by only 6% of discs hosted in halos with 
M = 2 x 1O 7 M0. The results obtained assuming md = 0.1 
for M = 10 7 M Q (solid line) and M = 0.7 x 1O 7 M (dotted 
line) are shown in Fig. [5] In both cases the possibility of 
fragmentation is included, but it only affects the lowest spin 
parameter halos, for which the probability is small to start 
with. The two panels refer again to two different choices of 
Q c . As expected, in this case, where a larger fraction of the 



halo mass collapses to the disc, the fraction of gravitation- 
ally unstable discs is larger, resulting in a larger probability 
of hosting a large mass concentration at the center. 

4.2 The case of molecular hydrogen cooling 

We now consider the case where molecular hydrogen can 
form and the disc therefore cools down to T gas w 500 K. 
In this case, discs forming out of haloes with T v i r w 10 4 
K will most likely fragment. However, an important point 
to note is that both the conditions for fragmentation, and 
the fraction of halo mass that will accumulate at the cen- 
ter only depend on the ratio T gas /T v i r as shown in eqn. 
|Q. Therefore, the very same central mass accumulation m a 
that can be obtained for a 4000 K disc can be obtained 
for a disc with a temperature ten times lower, provided 
that the halo virial temperature is reduced by the same 
amount (i.e., from T vir » 10 4 to T vir » 10 3 ). The typical 
redshift for the collapse of such haloes is z ~ 20, higher 
than the z ~ 10 assumed above. Since the halo mass scales 
as M oc T 3 / 2 (l + z)~ 3//2 , the same results shown in Figs. 
and|21would still hold for a T gas w 500 K disc, provided that 
the mass scale (of both the halo and of the central mass) 
is reduced by a factor ~ 100. The typical mass of a cen- 
tral concentration in this case is of the order of 10 3 Mq. Our 
proposed mechanism might therefore also be relevant to the 
conditions just prior to the formation of the first stars. 

5 TIME-DEPENDENT MODELS 

In this section, we develop some time-dependent models to 
confirm the results obtained with the simple estimates in 
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Figure 2. Cumulative distribution of the mass available for the formation of a seed black hole. The parent haloes are at a redshift z = 10 
and we assume here that raj = 0.1. The line styles refer to haloes with mass M = 10 7 Mq (solid) and M = 0.6 X 1O 7 M0 (dotted). 



previous sections. The evolution of the surface density E in 
a viscous disc with a flat rotation curve is given by: 

r)Y 1 r) 2 

M=R8R?^ + ^ < 14 > 

where v is the viscosity and Ej n is a source term representing 
the inflow from the hal o. 

It has been shown (iGammiefeOOlllLodato fc Ricefeool) 
that the transport associated with the development of such 
gravitational instabilities can be well described in terms of 
an effective viscosity, i.e. within the framework of the stan- 
dard a description of accretion discs. It is then convenient 
to adopt the simple prescription fo r viscosity induced b y 
gravitational instability described bv lLin fc Pringlel £l990): 

v = a g c B H = a e — — , (15) 

7TG2-/ 

where a g , the stress provided by gravitational instabilities, 
is modeled as: 

^(f-l). (16) 

In the previous equation, Q c is a critical value of Q, below 
which gravitational instabilities are able to provide a stress 
in the disc, and 77 is a parameter that essentially determines 
how far from the critical value of Q the disc has to be in 
order to deliver a given stress. In practice, the general be- 
haviour of such discs is that they attain a given value of 
Q, such that the stress provided through a g is the value re- 
quired to pass on the incoming flux Ei n . Modifying the value 
of 7] merely changes the equilibrium value of Q at which this 
stress is provided. Effectively, apart from small changes to 
the normalization of E, the specific value of r\ does not af- 
fect our results, therefore confirming that the final state of 
the disc and the central accumulation of matter does not 



depend on the specific viscosity mechanism. We have gener- 
ally chosen 77 = 0.1, but have also explored other values and 
do not notice any appreciable difference. 

The mass input is modeled in such a way that the accre- 
tion rate from the baryons contained in the halo onto the disc 
Af h = 2-7T J £i n i?d-R is given by eqn. J5J for t < tg = 7-200 /Vh 
and vanishes for t > tg . In this way, the total mass accreted 
onto the disc is a fraction ma of the halo mass. The ra- 
dial dependence of Ei n is determined under the constraint 
that the total accreted angular momentum is a fraction ja 
of the halo angular momentum, with the assumption that 
ja = ma- This constraint, however, still leaves some freedom 
in the choice of E. For example, the same result is obtained 
in the two extreme cases where all the mass is added to 
the disc as a 5-function centered at Ro = \/2Ar2oo, or the 
case where Ei n has an exponential dependence on radius 
oc exp(— R/Ro), with Ro = \r2oo/V%- 

The results of two of these time-dependent calculations 
arre shown in figs. [3] and 2]. Both of them refer to the follow- 
ing choice of parameters: A = 0.02, ma — 0.05, c s /Vh = 0.7, 
Q c = 3. Fig. prefers to the case where the matter is added 
as a 5- function, while Fig. prefers to the case where an ex- 
ponential profile of Ei n is adopted. The left panels show the 
surface density in the disc at t = 0.7, 0.8, Its and at the end 
of the simulation, when all the matter has collapsed to the 
disc and the disc becomes stable and is not accreting. The 
middle panels show the profiles of a g at t = 0.7, 0.8, Its. 
The right panels show the time evolution of the disc mass 
(solid line) and of the total mass accreted from the halo 
(dot-dashed line) . The fraction of halo mass accreted in the 
center at the end of the simulation was m a w 0.00688 for the 
5-function case and m a w 0.0073 for the exponential case, 
in very good agreement with the simple estimates obtained 
from eqn. {Jjj. 
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Figure 3. Results of a simulation where the mass is added to the disc with a 5-function, with A = 0.02, = 0.05, c s /Vh = 0.7 
and Q c = 3. Left Panel: Surface density profiles (in units of 160Mq/pc 2 , assuming 1 gas 

= 4000 K and z = 10) at t = 0.7,0.8, 1% 

(solid, short-dashed and long-dashed lines, respectively) and at the end of the simulation (dotted line). Middle Panel: Profiles of a g at 
t = 0.7,0.8, Itff (solid, short-dashed and long-dashed lines, respectively). Right Panel: Time evolution of the disc mass (solid line) and 
of the mass collapsed onto the disc (dot-dashed line). 
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Figure 4. Results of a simulation where the mass is added to the disc with an exponential profile, with A = 0.02, raj = 0.05, c s /Vh = 0.7 
and Q c = 3. Left Panel: Surface density profiles (in units of 16OM0/pc 2 , assuming T gas = 4000 K and z = 10) at t = 0.7,0.8,1% 
(solid, short-dashed and long-dashed lines, respectively) and at the end of the simulation (dotted line). Middle Panel: Profiles of o g at 
t = 0.7,0.8,1% (solid, short-dashed and long-dashed lines, respectively). Right Panel: Time evolution of the disc mass (solid line) and 
of the mass collapsed onto the disc (dot-dashed line). 



A comparison of the evolution of disc mass for two cases 
with different values of A = 0.02 (black line) and A = 0.01 
(dashed line) in shown in fig. |^| The other parameters were 
set as before and the mass input function was a ^-function. 
The accreted mass at the end of the simulation with A = 0.01 
is m a = 0.019, once again in very good agreement with 
eqn. Q. 

An interesting feature that appears from these time- 
dependent models, is the value of the viscosity parameter a g 
during the simulation. If the disc were to deliver a steady 
mass accretion rate of Mh, as given by eqn. ©, it would 
require an a e = md/2(Vh/c s ) 3 « 0.07 for the chosen param- 
eters (comparing eqn. JHJ with eqn. HOIl ). However, since the 
system is not in a steady state, the actual value of a g needed 
to provide the required mass flux is actually lower, peaking 
at roughly a g « 0.04. This is quite important to determine 
whether we expect the disc to fragment and form stars or 



not. As described earlier, we expect fragmentation whenever 
the stress needed to transport the mass flux through the disc 
is larger than a threshold value a c m 0.06. With the param- 
eters chosen above, that requires a w 0.07, we would then 
expect fragmentation, whereas in fact an evolving disc will 
always have a lower value effective of a and is therefore not 
expected to fragment. 

The above behaviour can be understood quantitatively 
as follows: if some matter SM is delivered at a radius Rq, 
only a fraction / = (1 — Rq/Ra) will flow into the inner 
disc, while the remaining 8M(Ro/Rd) will be carried out 
to compensate for the angular momentum lost by infalling 
matter. Of course, this will leave an unbalanced mass flux 
in the outer disc that will lead to an increase of R^. This 
however happens on a very long time-scale, especially since 
the viscosity in the outer disc becomes very small. Using 
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Figure 5. Comparison of the time evolution of the disc mass for 
two simulations with different values of the halo spin parameter: 
A = 0.02 (solid line) and A = 0.01 (dashed line). The dot-dashed 
line shows as in Figs. Inland HI the total mass collapsed onto the 
disc. 



eqn. ||3J to obtain the outer disc radius and assuming Ro 
V^Xr^oo, we obtain: 



/= 1- 



V2X 



r200 



Ra 



m d 



(17) 



where we have also assumed j'd = m^. For the case consid- 
ered in Fig. El we obtain / w 0.57 which then lowers the 
required stress from a ~ 0.07 to a fa 0.04, in perfect agree- 
ment with the results of the simulation. We have therefore 
demonstrated the result derived in Section |3] that the mass 
flow into the inner disc during the disc build up is actually 
smaller than its steady-state value by a factor /. 



6 LUMINOSITY AND OBSERVABILITY 

In this section we estimate the luminosity and potential ob- 
servability these pre-galactic discs. Here we only concentrate 
on the hot discs dominated by atomic hydrogen cooling. The 
dissipation rate per unit surface area determined by the ac- 
cretion process is: 



D{R) = vE(RCl') 



M 
2^ 



(18) 



where we have assumed a flat rotation curve at Vh- The total 
luminosity emitted by the disc is: 

/•flout 

Ldisc = 2tt / RD(R)dR 



MV h 2 



flout ^ 



flill 



m d -^-\n(R out /R in ), (19) 



where R m and 7? out are the inner and outer disc radius, re- 
spectively, using M » maVjf/G. The outer radius of the disc 
is given by i? ou t = Rd ~ 100 pc, for a gravitationally unsta- 
ble disc, with A = 0.05 at z — 10. The inner radius is harder 
to estimate. Our model will not be valid in the innermost 
parts of the disc, where the influence of a growing black hole 
will become important. Inside this radius the rotation curve 
will start to rise in an approximately Keplerian way, and the 
disc will therefore become hotter. This inner limiting radius 
can be estimated from: 



,, GMbh , _ 1n _ 2 / M BH \ /10km/sec\ 
-Rin = — — — « 4.5x10 ( , , ) I tt I pc.(20) 



\10 3 MqJ V li 



Substituting these values into eqn. 119L we find that the 
luminosity arising from the outer part of the disc is: 



4000L 



( 



o 



V10km/s 



(21) 



Note that, while in principle eqn. 12H would tell us that 
larger values of Vh would provide a much larger luminosity, 
this is actually unlikely, since for larger Vh the disc would ef- 
ficiently fragment into stars rather than accrete at a higher 
rate, as discussed above. We therefore argue that the ac- 
cretion luminosity in these discs is unlikely to exceed the 
estimate above. Since the gas cools mainly through atomic 
hydrogen line emission, most of this luminosity will be emit- 
ted in the Lyman series, and in particular in Lya. 

In fact, most of the luminosity of the system will be re- 
leased in the inner disc, as the matter falls to the bottom of 
the potential well to ultimately form and feed the growing 
black hole. A detailed description of this process is beyond 
the sco pe of the present paper and has been discussed else- 
where dVolonteri fc Reed 1200.3 iBegelman et al]l200rf) . In- 
deed, the accretion rate onto the black hole or onto the 
'quasi-star' described in such models is one of their input pa- 
rameters. This parameter can be computed self-consistently 
in a time-dependent way based on our calculation presented 
in section|S] The results are shown in Fig.|!|]and Fig.|S] The 
upper panel of Fig. |S] shows the time evolution of the mass 
of the growing central object Mbh- The middle panel shows 
the evolution of the mass accretion rate onto it (solid line) 
and the corresponding Eddington rate for a black hole of 
mass Mbh with an efficiency of matter-luminosity conver- 
sion set equal to 0.1. As can be seen the initial phases are 
characterized by a strongly super-E ddington accretion phase 
(as expected, and as predicted by IVolonteri fc Reesl 120051 : 
IBegelman et al.ll2006t) . At later times, however, as the mass 
of the central object grows to reach roughly its asymptotic 
value, the accretion rate decreases, becoming sub-Eddington 
when the black hole has a mass of roughly 8 x 1O 4 M0 in this 
case. The bottom plot shows the luminosity expected to be 
radiated from such a system, assuming that during the early 
super-Eddington phase the output luminosity is still limited 
by the Eddington value, while in the subsequent phase it 
is limited by the decrease in the mass supply. The peak of 
the luminosity is expected to occur just before the transi- 
tion from the super- to the sub-Eddington regime and it 
corresponds to a few times W 9 Lq. 
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Figure 6. Upper Panel: mass of the central concentration as a 
function of time for the case shown in Fig. [3] Middle Panel: M in 
the inner disc (solid line) and the Eddington rate for a black hole 
of mass Mbh (dotted line). Bottom Panel: expected luminosity as 
a function of time, assuming that it is limited by the Eddington 
value at early times and by mass supply at late times. 

6.1 Prospects for observing these pre-galactic 
discs 

Given the low intrinsic luminosity of these pre-galactic discs 
(~ 4000 1/© ), the prospects for direct detection are not 
promising. However, including the accretion luminosity of 
the black holes that are being assembled in the centers of 
these discs, they are likely to be detectable by the James 
Webb Space Telescope (JWST). The best strategy for de- 
tecting these objects is to exploit their lensing by foreground 
clusters. Using massive, foreground cluster-lenses as gravita- 
tional telescopes, we can expect to detect highly magnified 
pre-galactic discs. Therefore, a survey through cluster lenses 
by JWST is likely to be the most viable detection strategy. 



7 DISCUSSION AND CONCLUSIONS 

In this paper, we demonstrate that massive central concen- 
trations can form naturally from pre-galactic discs that as- 
semble in dark matter halos at high redshift. The masses of 
these concentrations depend on key properties of the host 
halo, mass, spin and gas cooling in the halo. In particular, 
low spin halos and massive halos are most efficient in con- 
centrating gas in their centers. However, not all halos will 
be able to accrete gas into their centers. Taking into ac- 
count the possibility of fragmentation furthers restricts the 
formation of eventual seed black holes. Using simple stabil- 
ity criteria, we predict that fragmentation and subsequent 
star formation are the fate for gas in some fraction of haloes. 
We derive 3 interesting regimes that are determined by the 
ratio of T v ir/T gas , (i) if this ratio is greater than 2.9, haloes 



will form discs that will fragment and form stars and not 
directly form black holes; (ii) if this ratio lies between 1.8 - 
2.9, haloes with low spin will lead to central mass concen- 
trations; however, the lowest spin cases might be affected by 
fragmentation; (iii) if this ratio is lower than 1.8 the haloes 
will not produce fragmenting discs but will successfully ac- 
crete gas in their centers. Calculating the accretion rates 
for the fate of the accumulated gas we find super- Eddington 
rates leading to eventual black holes of masses up to 10 Mq. 

The model pro posed here has two i mport ant predeces- 
so rs in the work of lEisenstein fc Loebl jl995l) and in that 
of iKoushiappas et alJ B2004I) . but also presents significant 
differences with respect to both of them. Primordial gas 
acquires angular momentum t hrough tidal interaction with 
the surrounding jPee bles 1969j). This is generally measured 
through the parameter A, the distribution of which can 
be d etermined from numerical simulations iWarren et alJ 
1992). Thus, the centrifugal barrier is the main obstacle to 
the formation of any compact object in the center of pri- 
mordial galaxies. It is therefore not surprising that all mod- 
els (including our own) predict that the most favourable 
sites for black hole formation are haloes with low spin. It is 
also not surprising that any black hole formation model has 
to deal with the problem of angular m omentum removal. 
This issue was tackled in the work of lEisenstein fc Loebl 
il995tl . who calculated the viscous time-scale for primor- 
dial thin discs and assumed that black hole formation would 
only occur when the discs are so compact that their vis- 
cous time-scale is less than the typical time-scale for star 
formation. They found that only very rare haloes, with A 
smaller by at least a factor 30 with respect to the average 
spin, would be able to produce black holes. This is sub- 
stantially different from what we find here. Indeed, we find 
that even haloes with A as large as 0.02-0.03 do produce 
large central mass concentrations. There are two main as- 
pects iWith_jejp_ect_towhich_our model differs substantially 
from lEisenstein fc Loebl lll995|). The first one is the na ture 
of the viscosity mechanism. lEisenstein fc Loebl il995T) ex- 
plicitly neglect the contribution of gravitational instabili- 
ties and assume that viscosity is driven by "turbulence". 
In standard accretion di scs, turbulence is thoug ht to arise 
from MHD instabilities (iBalbus fc Hawlevl ll998). However, 
for such primordial discs, where the primordial magnetic 
field is very weak and the gas is predominantly atomic or 
molecular, it is unlikely that MHD instabilities would be 
dynamically imp ortant. The second , and m ost important, 
difference is that lEisenstein fc Loebl (^395) considered rel- 
atively thin discs, with H/R m 0.03, whereas our discs are 
substantially thicker. Since the viscous timescale is propor- 
tional to (H/ R)~ 2 , this means that the viscous timescale as 
estimated bv lEisenstein fc Loebl <ll99fj) is much larger than 
in our case. This in turn, is what leads to their more pes- 
simistic estimate of the rarity of black hole forming haloes. 
Fi nally, a s ignificant di fferen ce between our model and that 
of lEisenstein fc Loebl lll995t) is that our model naturally 
leads to a robust determination of the seed black hole mass 
function. 

IKoushiappas et al.l <l2004l) propose that black hole seeds 
with masses of the order of 1O 5 M (i.e. very similar to those 
obtained here) can form out of low angular momentum ma- 
terial in massive haloes. They assume that the lowest an- 
gular momentum material within the halo forms a compact 
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disc which is gravitationally unstable and accretes onto the 
center due to the effect of an effective viscosity driven by the 
instability, in a way similar to what we propose here. How- 
ever, they do not consider self-consistently the evolution of 
the surface density profile induced by the assumed viscosity 
mechanism. In their picture, the disc surface density sim- 
ply grows linearly with time (at the same rate, independent 
of the radius) and at any given time it reflects the original 
angular momentum distribution of the gaseous component 
of the halo. On the contrary, what determines the surface 
density profile is actually the viscosity mechanism (which, 
being related to gravitational instabilities, in turn only de- 
pends on Q). As we have shown in Section |^| (Figs. [3] and 
2J, where we numerically follow the evolution of such discs, 
the final E profile (and the total accreted mass m a ) is the 
same whether we add mass to the disc with an exponential 
profile or with a 5-function, representing two extreme cases 
in the original angular momentum distribution within the 
halo. As a consequence of this, the esti mates of black hole 
masse s and their distribution given by iKoushiappas et"al] 
(2004) artificially depend on the initial angular momentum 
distribution and also on the assumed viscosity law. On the 
other hand, we have also clearly shown that, as long as vis- 
cosity is driven by gravitational instabilities, the black hole 
mass distribution is independent of viscosity. 

In the present paper, for illustrative purposes, we have 
taken the simplifying assumption that these primordial discs 
live in the potential well of a simple isothermal sphere, with 
a given circular veloci ty Vh- For a more rea listic density pro- 
file, the NFW profile dNavarro et aljll997t) . in its innermost 
parts, the disc mass might dominate the potential well and 
become radially unstable, giving rise to bar-like instabili- 
ties. This, however, will leave our results and conclusions 
unaffected. Firstly, the instability crite rion for such global 
instabilities dChristodoulou et alj ll995) is rather similar to 
our adopted criterion based on Q, and marginal stability 
is found to occur for equivalent values of Q c w 2 — 3, as 
adopted here. Secondly, since the disc mass is only a small 
fraction of the halo mass, only the innermost parts of the 
disc will be subject to bar-like instabilities, at radii much 
smaller than the typical disc radius . The development of 
such instabilities might in fact enhance the accretion rate in 
the inner disc, but will not change the estimates of the total 
accreted mass obtained here. The important and interest- 
ing consequence of our model is that black hole masses of 
10 9 Mq powering the luminous SDSS quasars at z — 6 can 
form comfortably within the available time of 1 Gyr in the 
concordance cosmology from our seed masses of 1O 5 M at 
z ~ 10. 
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APPENDIX: THE STRUCTURE AND 
STABILITY OF SELF-GRAVITATING 
ZERO-METALLICITY ACCRETION DISCS 

We consider the structure and stability of self-gravitating 
discs embedded in the external potential generated by a dark 
matter halo. If the disc is cold enough (or massive enough) it 
can become self-gravitating, i .e. it will develop gravitational 
instabilities. It is well known jGammidl200ll:lLodato fc Ricel 
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Figure 7. Structure of self-regulated accretion discs with cooling dominated by atomic hydrogen. The plots show (clockwise from top 
left): the temperature T gas , the surface density £, the optical depth r and the viscosity coefficient a . The lines refer to the following 
pairs of (M, V h ): solid line (10 -2 , 10), dotted line (10~\ 10), dashed line (10 -3 , 10), dot-dashed line (10~ 2 , 100), where M is in units of 
Mq/yi and Vh in km/sec. 



2004, 2005) that once gravitational instabilities set in, the 
disc will rapidly evolve into a quasi-steady state of marginal 
gravitational instability, characterized by a constant profile 
of the Toomre parameter Q. 

We start by noting that the requirement that Q be con- 
stant, for our flat rotation curve disc, i mplies that if the disc 
is isothermal (as it is expected to be, see lOh fc Haimanl2002l 
and arguments below) then the surface density E oc i.e. 
the disc has the same density profile as a simple Mestel disc. 
This in principle, makes it very easy to include the contribu- 
tion of the disc itself to the radial gravitational field, since it 
also leads self-consistently to a flat rotation curve. However, 
we will not consider this contribution at present. The steady 
state structure of such Mestel radially self- gravitating accre- 
tion d iscs has been explored extensively bv lBertin fc Lodatd 
( 1999), who also include the modifications to the disc struc- 
ture expected to take place very close to the black hole, 
where its gravitational field starts to become important. 

Actually, it is possible to solve for the detailed disc 



structure by relaxing the assumption of isothermality of the 
disc, and by computing the temperature by requiring that 
radiative cooling (here assumed to be optically thin) is bal- 
anced by the heating provided by the accretion mechanism. 
This energy balance equation gives us a rela tionship between 
a and the cooling time t coo i JPringlelll98l) : 



1 



7(7 — 1) tcool ' 



(22) 



where 7 is the ratio of specific heats, td yn = R/Vh is the 



dynamical timescale, and t c 



nkT gse /n A(T gas ), where 



n is the number density, r ga s is the gas temperature and A 
is the relevant cooling function. 

Once the cooling function is provided, it is then possible 
to solve the basic disc equations describing self-regulation 
(eqn. (Q), angular momentum conservation (eqn. IIIOH ) and 
thermal equilibrium (eqn. 1221 ) for the three unknowns 
E(-R), Tg as (7?) and a(R), with the input parameters Q, M 
and Vh. 
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Figure 8. Structure of self-regulated accretion discs with cooling dominated by molecular hydrogen. The plots show (clockwise from 
top left): the temperature T gas , the surface density £, and the viscosity coefficient a. The lines refer to the following pairs of (M,V h ): 
solid line (10~ 2 , 10), dotted line (10 _1 , 10), dashed line (10~ 3 , 10), dot-dashed line (10~ 2 , 100), where M is in units of M Q /yr and V h in 
km/sec. 



We consider here the case where the disc is made of pri- 
mordial gas, with no metal enrichment. The cooling process 
is then going to be provided by hydrogen. We consider two 
cases: (a) when molecular hydrogen formation is suppressed, 
so that cooling is dominated by atomic hydrogen, and (b) 
the case where molecular hydrogen is present and therefore 
dominates the cooling function. 

7.1 Cooling processes: Atomic hydrogen cooling 

We consider here th e coo l ing fu nction for pure atomic hy- 
drogen provided by ICenl jl992l) (we have also consid ered 
a different cooling function, after iKatz fc Gunnlll99ll . and 
found no significant difference). The results of the calcula- 
tions are shown in Fig.|7|where we plot the resulting profiles 
of E, T gas and a for different values of the input parameters. 
Here, we have assumed Q = 1 and varied M and Vh. As can 
be seen, due to the steepness of the cooling function, the 
disc turns out to be almost isothermal, with the equilibrium 
temperature only very weakly dependent both on radius and 
especially on input parameters, with T gas ~ 4000 K (note 
that changing the mass accretion rate by two orders of mag- 
nitude, only modifies T sas by roughly 14%). In Fig. Q we 
also plot the optical depth of the disc r = Ekr, where the 
Rosseland mean opa cities kr for a primordia l composition 
has been taken from iMaver fc Duschil j2005l) . We see that 
the assumption of optically thin cooling is indeed valid over 
the whole disc. 



extremely low. We can see that the equilibrium tempera- 
ture is going to be much lower than in the case of atomic 
hydrogen cooling, as expected. 

Since the cooling function is less strongly dependent on 
temperature for molecular hydrogen cooling, the "thermo- 
stat" that keeps the temperature of the disc constant in the 
case of atomic cooling is less efficient, and the equilibrium 
temperature in this case depends more strongly on the input 
parameters. The equilibrium temperatures range from 200 
K to 700 K. The disc is still close to isothermal, with the 
temperature rising in the inner disc. 



7.2 Cooling processes: Molecular hydrogen 
cooling 

We consider next the situation wherein molecular hydro- 
gen formation is not suppressed and we assu me that it can 
form up to a molecular fraction xh 2 = 10 -3 llOh fc Haimanl 
2002) . The cooling f unctio n for molecular hydrogen is taken 
from iGalli fc Pallal (Il99§l) . Fig. QU shows E, T gas and a for 
this case. We do not plot the optical depth for molecular 
hydrogen cooling since the available opacity tables for pri- 
mordial composition do not cover the very low temperature 
range predicted. At any rate, the opacity is expected to be 



